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ABSTRACT 

We investigate the correlation between far-infrared (FIR) and radio luminosities in dis- 
tant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large 
Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer Cam- 
erA (LAB OCA), the Very Large Array (VLA) and the Giant Metre- wave Radio Telescope 
(GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 
25 0-/im- selected galaxies, we re-measure the 70-870-/im flux densities at the positions of 
their most likely 24-/nm counterparts, which have a median [interquartile] redshift of 0.74 
[0.25, 1.57]. From these, we determine the monochromatic flux density ratio, #250 (= logio 
[^250^111 /£i,400MHz]), and the bolometric equivalent, qi^. At z w 0.6, where our 250-/im fil- 
ter probes rest- frame 160-/im emission, we find no evolution relative to qieo for local galax- 
ies. We also stack the FIR and submm images at the positions of 24-fim- and radio- selected 
galaxies. The difference between seen for 250-/im- and radio- selected galaxies suggests 
star formation provides most of the IR luminosity in <100-/iJy radio galaxies, but rather 
less for those in the mJy regime. For the 24-/im sample, the radio spectral index is constant 
across < z < 3, but #ir exhibits tentative evidence of a steady decline such that qm 
oc (1 + ^)-° 15±0 03 _ significant evolution, spanning the epoch of galaxy formation, with 
major implications for techniques that rely on the FIR/radio correlation. We compare with 
model predictions and speculate that we may be seeing the increase in radio activity that gives 
rise to the radio background. 

Key words: galaxies: evolution — infrared: galaxies — radio continuum: galaxies. 



© 0000 RAS 



Ivison et al. 



1 INTRODUCTION 

The correlation between FIR and radio luminosities (e.g . 
van der Kruitlll97ll : IPickev & Salpeterlll984l : IdeJong et al.l ll985: 



Helou etalJll985h is believed to be due to a common link with 



massive, dust-enshrouded stars. During their brief lives these stars 
warm the dusty molecular clouds in which they were born, expel 
dust in their winds, then create radio-luminous supernova remnants 
(which are believed to accelerate cosmic- ray electrons), pe rhaps 
generating more dust during this finale (e.g. lDunne et al.l l2009b). 

The empirical FIR/radio correlation is regularly exploited 
in a variety of ways - to calibrate th e relationship between ra - 
dio luminosity and star-formation rate (Condon 1992; Bell 2003), 
for e xample, or to estim ate the distance to luminous starbursts 
(e.g. Carilli & Yun 1999), or their dust temperatures (Td - e.g. 
IChapman et al .112005k or to define samples of ra dio-excess active 
galactic nuclei (AGN - e.g. iDonlev et al ] |2005l) . The correlation 
has thus become a cornerstone of modern astronomy. It is clearly 
important to know whether the correlation breaks down at extreme 
luminosities, or varies with redshift, perhaps due to variations in 
magnetic field strength (e.g. lBerneTe t al. 2008), IR photon density, 
initia l mass function (IMF), dust composition or cosmic-ray flux 
fe.g. lRengaraianll2005l) . 

Recent efforts have concentrated on determining the re- 
lationship between radio and 24-/xm flux densities, ^24 = 
log 10 (524/xm/£i,4ooMHz), taking advantage of Spitzer's sensitiv- 
ity to h ot dust emission from very distant st ar-forming galax- 
ies (e.g. lAppleton et all |2004 llbar et al.l l2008h . A galaxy's 24- 
/nm luminosity is not a particularly reliable tracer of its FIR lu- 
minosity, however, being prone to uncertain contamination by 
AGN continuum emission and by spectral featur es due to silicates 
and poly cyclic aromatic hydrocarbons (PAHs - IPope etalJbOOd : 
iDesai et al.ll2007l). 

lAppleton et al.l noted that 'Ideally, it would be better to mea- 
sure a bolometric q, but insufficient data are available at longer 
waveleng ths to do this relia bly'. Here, we exploit new data from 
BLAST toevlinetalJl2009h and t he LABOCA ECDFS Submm 
Survey (LESS - Iweiss et al]|2009h to explore the correlation be- 
tween bolometric IR luminosity and if -corrected radio luminosity. 
We do this separately for FIR-, 24-/xm- and radio- selected galaxies. 
The latter catagories are explored in £14.41 stacking at the positions 
of 24-/xm and radio emitters. 

Before we describe our analysis, the question of how best to 
explore the FIR/radio correlation for FIR- selected galaxies merits 
discussion. A thorough treatment of the IR and radio properties 
of FIR- selected galaxies must be able to deal with the potentially 
FIR-loud and radio-weak fractions of the sample. If the available 
FIR data were well matched in resolution to typical optical/IR or 
radio imaging then comparing FIR and radio properties would be 
trivial. As things stand, however, the BLAST 250-/xm beam cov- 
ers ^500 x more area than a typical VLA synthesised beam and 
matching FIR and radio sources is non-trivial ( 33. 11 ). 

We must consider how to deal with cases where counterparts 
to our FIR emitters cannot be found. This leaves us with no way 
to determine their positions unambigiously, or even their FIR flux 
densities since they are often blended, or their radio flux densities 
(or appropriate limits), or their redshifts. We can not begin the anal- 
ysis by pinpointing the FIR emitters in the radio waveband since 
this intr oduces a strong bias in favour of the most luminous radio 
emitters. iMarsden et alJd2009|) showed that the 24 -um galaxy pop- 
ulation detected by Spitzer ( Magnelli et al.ll2009h can account for 
all of the FIR background, as measured by the Cosmic Background 



Explorer jPuget et al] 1 19961 : iFixsen et all 1 19981) . We have already 
noted that the correlation between 24-/xm and FIR luminosities is 
poor, but every FIR source should have a 24-/im counterpart and 
our approach to pinning down the positions and radio properties of 
the p arent sample i s to adopt the 24- /xm identifications proposed 
by iDve et al] |2009) or, where iDve et al] opts for a radio identifi- 
cation, our own 24-/xm identifications. Armed with these positions, 
we determine the radio properties and redshifts of the appropriate 
counterparts, including appropriate limits for any that lack radio 
emission. 

Our paper is laid out as follows. In ^2] we introduce a 250- 
/xm-selected sample of galaxies from BLAST and the Spitzer, 
LABOCA, VLA and GMRT data used to determine their IR and 
radio spectral energy distributions (SEDs). In 33.1l we cross-match 
our FIR sample with radio emitters, demonstrating how difficult 
this will be for deep surveys with Herschel. We cross-match with 
the most likely 24-/nm counterparts ( 33. 2K enabling us to explore 
the FIR/radio correlation for our 250-/im- selected galaxies in $4] 
unbiased by radio selection. We also investigate the correlation 
for 24-/xm- and radio- s elected samples (34 .41 and FOK via the 
stacking technique (e.g. Ilvison et all 12007 al : iDunne et al . 2009a; 
Mar sdenetal] [2009), looking at the effect of K corrections using 
several SED templates. We state our conclusions in ^5] 

Throughout this paper we assume a Universe with Q m = 
0.27, 0\ = 0.73 and H Q = 71 kms" 1 Mpc" 1 dSoergel et al] 
l2007h . 



2 OBSERVATIONS 

Our primary dataset, known as 'BLAST GOODS South - Deep' 
(wher e GOODS is the Gre at Observatories Origins Deep Sur- 
vey - Dickinson et al. 2003) was taken during BLAST'S Antarc- 
tic flight in 2006 dDevlin et al.ll2009h and comprises map£] cov- 
ering the ECDFS at 250, 350 and 500 /xm. We combine them 
with panchromatic data, including the deepest available 70-160- 
pm imaging from Spitzer - from the Far-Infrared Deep Extra- 
galactic Legacy survey (FIDEL), with 870-/xm LABOCA imaging 
from the 12-m Ata cama Pathfinder Experiment (APEX) telescope 
(Weiss et al. 2009), with high-r esolution 1 ,400-MHz imaging from 
the VLA dMiller et al.ll2008l : lBiggs et ai1l2009h and with new 610- 
MHz imaging from the GMRT. This latter dataset allows us to in- 
correct the 1,400-MHz radio luminosity accurately out to z ~ 1.3, 
and beyond if the spectra are intrinsic power laws. Previous work 
assumes radio emitters share a common spectral index, a & —0.8 
(wher e S oc v a ) - a poor approximation, as shown by llbar etaD 
d2009l) . Together, these data are capable - prior to operations with 
Herschel - of determining the bolometric IR luminosities of distant 
galaxies. 

In order to take full advantage of the best available multi- 
wavelength coverage, we limit the areal coverage of this study to 
that of the deep VLA imaging described in 32.41 a region of ra- 
dius 17.2arcmin centred on 002000 = 03 h 32 m 28 s .3,(5j2ooo = 
— 27°48 / 30 // (see Fig. [TJ- This area coincides with the deepest of 
the BLAST surveys. 



1 Available at http : / /blastexperiment . info 
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Figure 1. BLAST 250-/im image of the central region of ECDFS, with 
>3-cr sources circled, a includes instrumental and confusi on noise com- 
bined in quadrature, where cr confusion /cr instrumental ~ 2 dMarsden et all 
2009). For the deepest imaging planned for Herschel, this ratio will be &25. 
Similarly significant sources at 350 and 500 /im are labelled with x and 
+ symbols, respectively. The region covered by the deepest VLA imaging, 
with relatively low levels of bandwidth smearing, is within the 17.2-arcmin- 
radius circle. 



2.1 BLAST 250-, 350- and 500-/im measurements 

The 1.8-m Balloon-borne Large Aperture Telescope, BLAST, is 
a forerun ner of the Spectra l and Photometric Imaging Receiver 
(SPIRE - iGriffin et al.ll2Q09h on the 3.5-m Herschel Space Obser- 
vatory. Data here are from the successful 2006 flight, lasting 270 hr, 
of which 90 hr was spent surveying a field centred on the ECDFS. 

In the part of this survey used here, the mean instrumental 
[confusion] noise levels are 11 [21], 9 [17] and 6 [15] mJybeam -1 
at 250, 350 and 500 pm, respectively, where the beams had FWHM 
of 36, 42 and 60arcsec ( Devlin et al J 120091 : iMarsden et alJl2009h 
and the noise values were determined in beam- smoothed maps. 

As our parent catalogue we adopt the samp le of 250-/xm 
BLAST sources described by Devlin et al ] (l2009h . cutting at a 
signal-to-noise ratio (SNR) of three, where the uncertainty in flux 
comprised the quadra ture sum of the instrumental noise (used 
to define SNR in the iDevlin et"aD catalogue) and the confusion 
noise. The reliability of our sample is thus s imilar to those of 
the first submm ga l axy (SMG) catalog ues (e.g. ISmail et aljfl997l : 
iHughes et al1ll998l : iBarger et al.lll998l) . Thus defined, our parent 
catalogue comprises 22 sources, listed in Table Q] and marked in 

Fig. in 

For the FIR measurements described in the sections that fol- 
low, the BLAST image in each of the three filters has been con- 
volved with the appropria te noise-weighted point spread function 
(PSF - iTruch et al.l l2009h . This is equivalent to calculating the 
maximum- likelihood point- source flux density by which the PSF 
would have been scaled to fit an isolated point source centred in 
each pixel. Flux densities in each band are taken from these con- 
volved maps at the appropriate positions. 

Since the differential 250 - /mi source counts f all rapidly with 
flux (dN/dS oc S~ 3 5 - see lPatanchon et alJl2009h . the effects of 
Eddington bias and source confusion have conspired to boost the 



fluxes in our parent catalogue. To determine the degree of boost- 
ing, we simulate these effects by i njecting a sourc e population 
consistent with the BLAST counts jPatanchon et al.ll2Q09h into a 
Gaussian noise field resembling the instrumental fluctuations in the 
BLAST GOODS South - Deep field and then measuring the flux 
densities of thos e sources at the input positions (see Appendix B of 
Eales et al . 2009). These simulations allow us to apply appropriate 
first-order corrections and to estimate the uncertainties associated 
with those corrections (see later in £14.21 ). The corrections are statis- 
tical in nature and can be tailored to individual sources only in the 
sense that a source with a given SNR is likely to have been boosted 
by factor, B (SNR), where B ranges from ~ 1.1-1.6 for the range 
of SNRs present in our sample. A full treatment of flux boosting in 
the BLAST data - including the influence of clustering and multi- 
band source selection - will be explored in future work. 



2.2 Spitzer 70- and 160-/im measurements 

All of the raw FIDEL (PL: Dickinson), GO (P.L: Frayer), and GTO 
(RL: Rieke) Spitzer data covering ECDF S were reduced and com- 
bined using the techniques developed by Fraver et al . (2004 We 
adopted the updated calibration factors, including corrections for 
colour and for the point- source-response function (PRF) given by 
lFraveretalJd2009h . 

The Spitzer flux density measurements are made at the po- 
sitions discussed later and the quoted uncertainties represent the 
r.m.s. noise at those positions, combined with an additional 10 (15) 
per cent uncertainty at 70 pm (160 pm) due to possible systematics 
affecting the calibration factors. 



2.3 LABOCA 870- pm measurements 

LESS covers around 900 arcmin 2 of ECDFS a t a wavelength of 
870 //m, with a beamwidth of 19.2 arc sec FWHM dGreve et alJ 2009: 
ICoppin et al1l2009l : IWeiss et al.ll2009h . to an average depth of a = 
1.2 mJy beam -1 as measured in a beam- smoothed map across the 
region in which our 250-/xm sample was selected. The 870-/im flux 
densities, at the positions discussed later, are measured in the same 
manner discussed in £12.1l and l22l 



2.4 VLA 1,400-MHz imaging 

Deep, high -resolution 1 , 400-M Hz imaging of the ECDFS was de- 
scribed bv lMilleretaH(l2008h and we use their image to identify 
radio counterparts for our parent 250-/xm- selected source catalogue 
(see g3j. 

We use the techniques described bv llbar et~al] d2009h to gen- 
erate a new, deeper 1,400-MHz catalogue and to correct for effects 
such as flux boosting at low SNRs. This catalogue was cut at > 4 a, 
where a was determined from the local background. 

In cases where the emission is found to be heavily resolved, 
e.g. for the large spiral associated with BLAST J033235— 275530, 
we measure the total radio flux density using TVS TAT within 
MVS. 
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Figure 2. Spitzer 24- /im greyscale images superimposed with red 1,400-MHz contours and blue 250-/im contours, centred on BLAST 250-/xm-selected 
sources which are labelled with their names. Each postage stamp image is 90arcsec x 90arcsec. Dashed black circles show our search areas (3 cr p0 s radius - 
see ^3). Broadly speaking, the sources are associated with either large spirals or clumps of faint, 24- fim galaxies - sometimes several of the former or as many 
as seven of the latter. The number of double sources is striking, especially those where the 250-/im emission peaks between the two. 



2.5 GMRT 610-MHz imaging 

To test the feasibility of deep radio observations in ECDFS, a small 
amount of new 610-MHz data were obtained using the GMRT0 
in 2008 November. During seven 9-h r sessions we obta ined data 
centred on the six positions chosen bv lMiller et alJ (2008) for their 
VLA survey f £12.4K recording 128 channels every 16 s in the lower 
and upper sidebands (602 and 618 MHz, respectively) in each of 
two polarisations. The integration time in each field was ^3hr, 

2 We thank the staff of the GMRT who made these observations possible. 
GMRT is run by the National Centre for Radio Astrophysics of the Tata 
Institute of Fundamental Research. 



made up of 30-min scans sandwiched between 5-min scans of 
the nearby calibrator, 0240—231, with 10-min scans of 3C48 and 
3C 147 for flux and bandpass calibration. 

Calibration followed standard recipes within AXVS 
(31DEC09). A bandpass table was generated after self-calibrating 
3C48 and 3C 147 in phase with a solution interval of 1 min. This 
was applied during subsequent loops of calibration and flagging 
(with UVFLG, TVFLG, SPFLG and FLGIT). The error- weighted 
flux densities of 0240-231 were found to be 5.55 ± 0.02 and 
5.31 ± 0.02 Jy at 602 and 618 MHz, respectively; the r.m.s. scatter 
of the flux density measurements over the seven days was 0.12 and 
0.10 Jy which suggests that the flux calibration should be accurate 
to better than 5 per cent. 



© 0000 RAS, MNRAS 000, 000-000 



Bolometric FIR/radio correlation 5 



The calibrated data were averaged down to yield 41 chan- 
nels in each sideband and the data from the seven days were con- 
catenated using DBCON. The resulting datasets were then self- 
calibrated in phase (SOLMODE=p!a), with a solution interval of 
2min, using CLEAN components from our own reduction of the 
VLA A-configuation data. This process yielded accurate phases; 
the astrometric reference frame was also tied to that of the 1,400- 
MHz image as a result. Subsequent imaging entailed the creation of 
a mosaic of 37 facets - to cover most of the primary beam in each 
of the six pointings - each facet with 512 2 pixels (1.5 2 -arcsec 2 per 
pixel). A further 10-20 bright sources outside these central regions, 
identified in heavily tapered maps, were also imaged for each point- 
ing. Subsequent self-calibration was carried out in phase alone, 
then in amplitude and phase (first with CPARM = 0, 1,0 then with 
CPARM = 0), with a solution interval of 2min, staggered by 1 min. 
The uv data were weighted using ROBUST = —0.5, UVRANGE = 
0.9, 100 kA, GUARD = — 1, UVTAPER = 30, 60 kA and UVBOX = 
10. 

After CLEAN components were subtracted from the uv data, 
more manual flagging was undertaken. Almost one third of the 
baselines were rejected, in total. CLEAN components were re- 
introduced (UVSUB, FACTOR=— 1), then the final six mosaics were 
convolved to a common beam size (6.5arcsec x 5.4arcsec, with 
the major axis at position angle 174°), then knitted together using 
FLATN. An appropriate correction was made for the shape of the 
primary beam, with data rejected at radii beyond the half-maximum 
level. The final image has a noise level of ^40/iJy beam -1 and 
covers the entire VLA image, with insignificant levels of bandwidth 
smearing. We exploit it here to determine the 610-MHz flux den- 
sities, and hence the spectral indices between 610 and 1,400 MHz, 
a 6io°°' °f me ra dio identifications described in ^3] These measure- 
ments are listed in Table Q] Our 610-MHz image is available on 
requesQ 



3 IDENTIFICATIONS 

In order to avoid introducing a strong radio-related bias, and despite 
the problems caused by the high surface density of 24-/xm emitters 
(e vident in Fig. Q, we use the 24- /am counterpar t identifica tions 
of iDve et all (120091) for our parent sample. Where IDve et al] does 
not present an association, or where they favour a radio identifica- 
tion over a 24-/xm source, we adopt our own 24-/xm identifications 
based on the probabilistic approach outlined in the next section. 

Before we discuss the 24-/xm identifications - and although 
we do not use them hereafter - we first describe the identification 
of radio counterparts to the parent FIR catalogue. Radio imaging 
remains the 'gold standard' for pinpointing submm and FIR galax- 
ies because of the tight flux correlation and the low surface density 
of radio emitters relative to optical/IR sources, and relative to the 
large submm beam sizes. It is important, therefore, that we identify 
any impending problems facing deep FIR surveys with Herschel, 
partic ularly in the SPIRE bands (250-500 /xm) where iDevlin et all 
(2009) predict confusion limits only 1.5 x lower than those suffered 
by BLAST. 

3.1 Radio identifications 

The positional uncertainty of a source found in a map is given by 

3 E-mail: rji@roe.ac.uk 
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Figure 3. Illustration of how the integrated value of P for identifications 
with P < 0.05 (so £P(P < 0.05)) for our sample varies with choice of 
search radius, where search radius is in multiples of cr p0 s (the r.m.s. offset 
between observed and actual coordinates in R.A. or Dec. - see equation[TJ. 
This is the solid line, and the left-hand axis. In a situation not dominated by 
the effects of confusion, the value at the chosen radius equates directly to 
the expected number of spurious identifications amongst our sample. The 
right-hand axis shows the total number of identifications in our sample, and 
the number that we would expect to miss by selecting too small a search 
radius. The dotted line shows the adopted search radius - in this case o (7 pos 
- where we have 12 formal identifications and ~0.24 was likely to have 
fallen outside our search radius. 



apos - 0.6 [SNR 2 - (2 7 + 4)]" 1/2 FWHM. (1) 

for a situation w here the counts ob ey a power law of the form 
N(> S) ex S' 1 dlvison et al.ll2007bh . The SNR term here is deter- 
mined optimally in a beam- smoothed map and includes any vari- 
ance due to confusion. The signal is uncorrected for flux boosting 
and the positional offset is not radial but is instead measured in 
Right Ascension or Declination. 

To associate radio emit ters with the FIR sour ces in our sample, 
we utilise the technique of Downes et al as implemented 

bv llvison et al.l d2002h and by most FIR/submm identification work 
thereafter. According to lDownes et all P is the probability of a ran- 
dom association and depends on the surface density of radio emit- 
ters and the radius within which one searches. The recent conven- 
tion has been to accept associations as secure identifications where 
P < 0.05. We assume a slope of 1.5 for the radio counts in the rel- 
evant flux density regime (^40-400 /xJy) and fine-tune the absolute 
level of the counts with Monte Carlo simulations. 

In choosing the search radius, we wish to maximise the num- 
ber of secure, unambiguous identifications and to minimise the 
number of real counterparts missed. Fig. [3] illustrates how the to- 
tal number of secure identifications and the integrated value of P 
for those sources [hereafter £P(P < 0.05)] varies with choice 
of search radius. £P(P < 0.05) equates directly to the expected 
number of spurious identifications in the sample. It is not intuitively 
obvious why there is a well-defined peak in SP(P < 0.05) when 
the total number of secure identifications is roughly constant with 
search radius: we see this because faint, real radio identifications 
compete (in terms of low P) with rare, bright sources - radio-loud 
AGN, entirely unrelated to the FIR source; as the search radius 
increases, correct identifications drop out of the integral (simply 
because their P rises above 0.05); the resulting reduction in se- 
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Table 1. Parent 250-/mi catalogue from BLAST (Devlin et al. 2 009h . 24- /j,m identifications, radio and X-ray properties, and redshifts. 



BLAST name c*25(Vm ^SO^m ^Iso^m a 24^m ^24^m S 2 4^m S 1)4 oomHz -^ghz ^610°° X-ray?* Spec Phc 

J2000 J2000 /mJy J2000 J2000 /mJy //xJy //xJy z z 



J033235 


-275530 


03:32:35.09 


— 27:55:31.0 


176.8 ± 10.9 


03:32:35.07 


—27:55:32.6 


4.7 


680 ± 34* 


565 ± 99 


±0.22 


Y 


0.038 


0.051 


J033229 


-274414 


03:32:29.74 


— 27:44:14.4 


156.7 ± 10.9 


03:32:29.87 


—27:44:24.2 


11.1 


1,558 ±55* 


1,670 ± 84 


—0.08 


Y 


0.077 


0.077 


J033250 


-273421 


03:32:50.01 


— 27:34:21.3 


159.3 ± 11.1 


03:32:50.41 


—27:34:20.3 


2.9 


470 ±20* 


706 ± 56 


—0.49 


N 


0.251 


0.251 


J033128 


-273916 


03:31:28.71 


— 27:39:16.1 


105.3 ± 11.1 


03:31:28.82 


— 27:39:16.7 


0.46 


34.8 ±7.5 


3a < 150 


—0.80 


N 






J033249 


-275842 


03:32:49.47 


-27:58:42.0 


101.2 ± 11.0 


03:32:49.34 


-27:58:44.7 


0.32 


216 ± 16 


296 ± 49 


-0.38 


N 




2.215 


J033124 


-275705 


03:31:24.05 


—27:57:05.5 


96.3 ± 11.1 


03:31:23.48 


—27:56:58.5 


1.3 


165 ± 16 


183 ± 52 


—0.12 


N 




2.732 


J033152 


-273931 


03:31:52.82 


— 27:39:31.5 


93.4 ±11.0 


03:31:52.07 


—27:39:26.6 


0.20 


965 ± 16 


765 ± 52 


±0.28 


Y 




2.296 


J033258 


-274322 


03:32:58.24 


— 27:43:22.3 


92.6 ± 11.1 


03:32:59.20 


—27:43:25.1 


0.53 


148 ± 17 


205 ± 43 


—0.39 


N 




1.160 


J033129 


-275722 


03:31:29.79 


-27:57:22.6 


91.6 ± 11.0 


03:31:29.92 


-27:57:22.4 


0.27 


144 ±16 


169 ±51 


-0.19 


N 


— 


1.571 


J033319 


-275423 


03:33:19.37 


-27:54:23.2 


90.7 ± 10.9 


03:33:18.89 


-27:54:33.8 


0.37 


60.2 ± 15.0 


3a < 150 


-0.80 


N 




0.488 


J033246 


-275744 


03:32:46.05 


-27:57:44.0 


90.0 ± 10.9 


03:32:45.88 


-27:57:44.8 


2.7 


337 ± 17* 


364 ± 89 


-0.09 


Y 


0.103 


0.105 


J033149 


-274335 


03:31:49.71 


-27:43:35.9 


86.6 ±11.0 


03:31:49.69 


-27:43:26.4 


0.95 


192 ± 16 


271 ±97 


-0.41 


N 


0.618 


0.603 


J033217 


-275905 


03:32:17.07 


-27:59:05.8 


85.3 ±11.0 


03:32:16.99 


-27:59:16.0 


0.37 


178 ± 17 


395 ± 45 


-0.96 


N 


0.126 


0.220 


J033318 


-274610 


03:33:18.13 


-27:46:10.3 


82.6 ± 10.9 


03:33:17.78 


-27:46:05.9 


0.43 


100 ± 14 


194 ±48 


-0.80 


N 




2.059 


J033216 


-280345 


03:32:16.08 


-28:03:45.0 


82.1 ± 11.0 


03:32:15.90 


-28:03:47.1 


0.77 


192 ± 20 


268 ± 49 


-0.40 


N 




0.452 


J033145 


-274635 


03:31:45.54 


-27:46:35.5 


80.2 ± 10.9 


03:31:45.36 


-27:46:40.2 


0.29 


42.2 ±7.3 


171 ±53 


-1.68 


N 






J033308 


-274805 


03:33:08.56 


-27:48:05.9 


79.7 ± 10.8 


03:33:09.72 


-27:48:01.3 


0.20 


433 ± 19 


496 ± 45 


-0.16 


N 


0.180 


0.190 


J033151 


-274431 


03:31:51.15 


-27:44:31.8 


74.0 ± 10.9 


03:31:51.09 


-27:44:37.1 


0.52 


96.3 ± 13 


3a < 150 


-0.53 


Y 




1.911 


J033130 


-275604 


03:31:30.07 


-27:56:04.5 


74.0 ± 10.9 


03:31:30.05 


-27:56:02.2 


0.90 


1,270 ±64* 


1,330 ±67 


-0.06 


Y 


0.677 


0.727 


J033218 


-275216 


03:32:18.03 


-27:52:16.8 


73.6 ± 10.9 


03:32:18.25 


-27:52:24.9 


0.21 


23.7 ±7.2 


3a < 150 


-0.80 


N 


0.739 


0.724 


J033240 


-280314 


03:32:40.11 


-28:03:14.7 


72.7 ± 11.0 


03:32:39.69 


-28:03:10.5 


0.34 


168 ±30* 


707 ± 180 


-1.73 


N 




0.956 


J033145 


-274205 


03:31:45.02 


-27:42:05.2 


72.8 ± 11.0 


03:31:44.46 


-27:42:12.0 


0.43 


122 ± 19 


181 ±45 


-0.47 


N 




1.056 



Notes: "^Flux densities before deboosting (see $\4.2\ ; thes e unce rtainties include only instrumental noise and should be added in quadrature with the 

confusion noise (c con f usion = 21.4 mJy - iMarsden et a l. 2009) to ar rive at realistic S NRs; *Total flux density, determined using tvstat. *24-/im position 
coincident - within 2 arcsec - of a Chandra source ( Lehmer et al. 2005; Luo et al. 2008). No further matches are made if the search radius is doubled. 



cure identifications is balanced by radio-loud AGN - contaminants; 
since they are very bright, their P tends to zero as soon as they are 
within the search radius; thus EP(P < 0.05) drops gently while 
the number of apparently 'secure' identifications is maintained. Al- 
lowing a search radius well beyond that of the peak P is therefore 
not advisable. 

We express the search radius in terms of multiples of <j pos to 
maintain an appropriate dependence on the SNR of the source for 
which a counterpart is sought. Based on Fig. [3] we search within 
a radius of 3 a pos , though values as low as 1.5 <j pos are defensible 
choices. This yields 12 statistically secure identifications with an 
expectation that ~0.24 lie outside the search area (calculated via 
the cumulative Rayleigh distribution function). 

At face value, the fraction of 250-/im sources with secure ra- 
dio i dentifications, 55 per cent, is similar to that found for SMG s 
fe.g. lSmail et al.l2000l : ivison et alJl200il2007bl : IPoDe et al.l2006h . 
which are expected to have higher luminosities and higher redshifts 
dChapman et al.|[2Q05h . However, close examination of Fig. |3 in 
which we have superimposed the radio and 250-/im data on a 24- 
fim greyscale image of each BLAST source, reveals that our abil- 
ity to securely associate radio counterparts with FIR emitters has 
been compromised. Broadly speaking, the BLAST sources are as- 
sociated with either large spirals or clumps of faint, radio and/or 
24-/xm galaxies - sometimes several of the former or as many as 
seven of the latter (e.g. BLAST J033 152-273931). At least half of 
the FIR sources appear to be blends of two or more faint galaxies - 
each often with faint radio identifications - i .e. confusion is a ma- 
jor problem dCondonlll974l ; lBlain et al .11 199 8). This is compounded 
by the large number of radio sources found by chance in the re- 
quired search area (a situation that worsens at 24 /xm, as mentioned 
earlier). Both issues are related to the relatively large size of the 
250-/xm beam (and the situation is worse at 350 and 500 /xm). 

As an alternative approach, we inspected the data shown in 
Fig. |3 by eye, rejecting seven sources with complicated 250-/im 
morphologies suggestive of multiple, blended emitters. Of the re- 
maining fifteen, nine have statistically secure radio identifications; 



of these, only three have unambiguous, single identifications. In to- 
tal, therefore, only 14 per cent of our parent sample have secure, 
unambiguous radio identifications, based on this more subjective 
approach. 

Thus, the purely statistical technique advocated previously for 
SMG identification work proves unhelpful - even misleading - in 
the confusion-dominated regime faced by BLAST and anticipated 
for the deepest survey work with Herschel (see Fig. Q] caption). 
Working through the S pitzer and Photodetec tor Array Camera and 
Spectrometer (PACS - IPoditsch etal.ll2008h bands to identify the 
dominant contributors to SPIRE flux densities may be more fruitful 
- but this has yet to be tested. This argues that a stacking approach 
may be the most enlightening - we turn to this later. 

3.2 24-/im identifications 

Using a different probabilistic approach to our own, iDve et al] 
(2009) listed the most likely 24-/xm counterparts for many 
of our 250-/xm- selected galaxies. To their list, we added 
BLAST J033229-274414, an obvious blend of t wo galax ies at 
z ~ 0.08 (Fig. |3, and in a further six cases where IDve et al.l listed 
no association we added the most likely of our own 24-/xm identi- 
fications. We list the 24-/xm positions in Table Q] together with the 
radio properties at those positions. The FIR and submm properties 
at the 24-/xm positions are given in Table |3 

As already noted, we adopt these 24-/xm counterparts to avoid 
introducing a strong radio-related bias, and despite the problems 
caused by the high surface density of 24-/xm emitters, exacerbated 
by FIR confusion. 



4 RESULTS 

4.1 Redshift distribution 

Of the most likely 24-/xm counterparts for our 250-jimi sam- 
ple, eight have spectroscopic redshifts (see ICroom et alJ l200ll : 
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Table 2. FIR-submm photometry at the positions of the 24-fim identifications, prior to deboosting, plus radio luminosities, temperatures and adopted redshifts. 



Name 




SVOjLim 


S'l60Aim 


<S250/xm 


S350/-i 






S870n 




T d (obs) 


Sir 


^l,400MHz 


Adopted 






/mJy 


/mJy 


/mJy 


/mJy 




/mJy 


/mJy 




/K 




/10~ 15 Wm -2 


AVHz -1 


z 


BLAST J033235 - 


275530 


93.9 ± 11.2 


351.1 ±65.0 


176.8 ± 10.9 


75.4 ± t 


.6 


29.8 ±6.0 


1.23 ± 1 


16 


25.2+0 


7 

9 


8.28 ± 1.79 


2.18 x 


10 21 


0.038 


BLAST J033229- 


274414 


136.6 ± 15.4 


368.0 ±65.8 


141.6 ± 10.8 


71.7 ±1 


.5 


36.8 ±6.0 


0.88 ± 1 


17 


27.2+5 


9 



10.5 ±2.8 


2.13 x 


10 22 


0.077 


BLAST J033250- 


273421 


30.1 ± 4.9 


203.5 ±43.1 


154.8 ± 11.0 


76.9 ± t 


.7 


40.3 ± 6.0 


1.82 ± 1 


22 


21.51° 


5 
9 


4.22 ± 1.24 


8.45 x 


10 22 


0.251 


BLASTJ033128- 


273916 


7.1 ± 1.9 


88.5 ± 25.2 


105.3 ± 11.1 


69.6 ± t 


.7 


39.8 ±6.3 


4.02 ± 1 


18 


"■it? 


6 
8 










BLASTJ033249- 


275842 


8.4 ± 2.1 


72.6 ±21.7 


101.2 ± 10.9 


66.4 ± t 


.6 


22.6 ±6.0 


1.75 ± 1 


13 


19.3+5 


7 




1.28 ±0.63 


6.38 x 


10 24 


2.215 


BLASTJ033124- 


275705 


8.7 ± 2.3 


100.7 ± 30.1 


83.2 ± 11.1 


47.2 ± i 


17 


33.3 ± 6.2 


2.99 ± 1 


20 


19.0+0 


7 

9 


1.40 ±0.70 


7.07 x 


10 24 


2.732 


BLASTJ033152- 


273931 


3.0 ± 1.4 


31.9 ± 14.4 


78.3 ± 11.0 


64.3 ± t 


.6 


53.1 ±6.0 


3.49 ± 1 


14 


15.9+1 


7 

8 


0.77 ± 0.42 


2.35 x 


10 25 


2.296 


BLASTJ033258- 


274322 


7.4 ± 2.0 


82.8 ± 24.8 


71.5 ± 11.1 


53.6 ± t 


.7 


34.3 ±6.0 


4.31 ± 1 


17 


18.3+0 


7 

9 


1.21 ±0.61 


9.48 x 


1Q 23 


1.160 


BLASTJ033129- 


275722 


3.6 ± 1.7 


49.0 ± 19.9 


91.6 ± 11.0 


54.4 ± i 


l.l 


46.3 ± 6.2 


5.11 ± 1 


14 


16.ll 1 , 


4 



0.93 ± 0.49 


1.78 x 


10 24 


1.571 


BLASTJ033319- 


275423 


3.9 ± 1.4 


36.0 ± 14.9 


84.2 ±11.0 


20 .4 ± i 


1.5 


13.2 ± 5.9 


1.12 ± 1 


13 


18.8+ 1 


■2 

4 


0.71 ±0.41 


5.24 x 


10 22 


0.488 


BLASTJ033246- 


275744 


46.0 ± 6.0 


147.9 ± 32.9 


90.0 ± 10.9 


46.7 ± I 


1.5 


18.7 ± 5.9 


1.30 ± 1 


14 


24.3+ \ 






3.03 ± 1.34 


8.49 x 


10 21 


0.103 


BLASTJ033149- 


274335 


15.7 ± 2.8 


100.5 ± 24.3 


81.5 ± 11.0 


46.6 ± i 


1.6 


30.9 ± 5.9 


2.46 ± 1 


16 


20.7^5 


8 



1.51 ±0.79 


2.75 x 


10 23 


0.618 


BLASTJ033217- 


275905 


39.5 ± 5.9 


110.4 ± 30.4 


80.5 ±11.0 


63.6 ± i 


1.5 


30.4 ± 6.0 


1.62 ± 1 


14 




1 
1 


2.00 ± 1.04 


7.17 x 


10 21 


0.126 


BLASTJ033318- 


274610 


5.6 ±1.7 


62.3 ± 20.8 


79.9 ± 10.8 


72.5 ± i 


1.6 


51.4 ±5.9 


3.84 ±1 


13 




8 
9 


1.00 ±0.56 


2.94 x 


10 24 


2.059 


BLASTJ033216- 


280345 


16.0 ±3.2 


85.7 ± 24.5 


82.1 ± 11.0 


39.2 ± i 


1.5 


-0.3 ±6.0 


2.04 ± 1 


30 


21.2+J 



1 


1.31 ±0.70 


1.31 x 


10 23 


0.452 


BLASTJ033145- 


274635 


-0.3 ±1.0 


12.0 ± 10.4 


75.2 ± 10.8 


11.0 db ^ 


1.5 


12.8 ±5.9 


0.30 ± 1 


14 


15.211 


1 
7 










BLASTJ033308- 


274805 


53.0 ±6.6 


97.9 ± 26.6 


66.1 ± 10.8 


37.0 ±i 


1.5 


31.7 ±5.9 


2.02 ± 1 


16 


25.9+1 


3 
5 


1.99 ± 1.04 


3.65 x 


10 22 


0.180 


BLASTJ033151- 


274431 


6.0 ±1.7 


77.5 ± 20.2 


74.0 ± 10.8 


63.8 ±i 


1.5 


37.4 ±5.9 


4.18 ±1 


16 


17.91° 


7 

9 


1.05 ±0.61 


2.13 x 


10 24 


1.911 


BLAST J033 OO- 


275604 


28.5 ±4.1 


83.2 ± 24.7 


74.0 ± 10.8 


20.9 ± i 


1.5 


44.9 ± 6.0 


3.05 ± 1 


15 


23.0+1 


1 

2 


1.61 ±0.91 


2.10 x 


10 24 


0.677 


BLAST J03321 8- 


275216 


5.3 ±1.4 


54.2 ± 17.8 


71.4 ±10.8 


52.6 ± i 


1.5 


46.8 ±5.9 


7.74 ± 1 


17 


16.ll 2 , 


3 
1 


0.88 ±0.52 


5.63 x 


10 22 


0.739 


BLASTJ033240- 


280314 


6.8 ±2.0 


50.7 ± 18.6 


72.7 ± 11.0 


61.8 ± J 


1.5 


30.6 ±6.1 


4.05 ± 1 


25 


18.211 


4 
4 


0.90 ± 0.54 


9.53 x 


10 23 


0.956 


BLASTJ033145- 


274205 


4.8 ±1.6 


45.0 ± 17.2 


56.0 ± 11.0 


27.0 ± i 


\.l 


28.4 ±6.2 


1.98 ±1 


14 


18.811 



5 


0.70 ± 0.45 


6.38 x 


1Q 23 


1.056 



Notes: 4QQMHz is the K-corrected radio luminosity using the measured a. Sir is deboosted; this process dominates the error budget (see ^2. H and ^4.2) . 

Spectroscopic redshifts are shown in italics. 
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Figure 4. The redshift distribution of the most likely 24-fim associations 
with our BLAST 250-/im sample. The dotted line shows the redshift dis- 
tribution of the 24-fim sources used for stacking in ^4.41 The subset of 
galaxies with spectroscopic redshifts are shown as a hatched area. Inset 
is a histogram of S'soo/xm/^so/xm, where we expect galaxies at higher 
redshifts to have higher ratios. Sources at redshifts higher than the median 
(z = 0.74) are shaded and there is a very clear tendency for these to have 
higher ratios, as we would expect. 



Le Fevre et all 120041 : ISzokolv et all 120041 : IVanzella et all [2008; 
Eales et al. 20 091) a nd another 1 1 have photometric redshifts 
dRaffertv et al.l l2009, in preparation). For every spectroscopic case, 
both types of estimate are available, and these agree reasonably 
well (Sz Y .m.s. ^ 0.05, with one outlier). We were unable to assign 
redshifts with confidence in two cases (see Tables Q]and[2]). 

The redshift distribution of 24-/im-identified, 250-/im- 
selected sources is shown in Fig. 01 The median is z = 0.74 with an 
interquartile range of z = 0.25-1.57, so slightly beyond the 70-/xm- 



Figure 5. Deboosted IR luminosity as a function of redshift. The dashed 
line is an arbitrarily scaled detection threshold for our adopted cosmology. 
Galaxies with spectroscopic and photometric redshifts are plotted as circles 
and stars, respectively. 



selected sample of ISvmeonidis et al ] d20Q9h which had a median z 
of 0.42. Inset in the same plot is a histogram of Ssoo/xm/S^so/xm- 
This ratio should increase with redshift as the SED peak shifts to 
longer wavelengths and we find a strong tendency for those above 
the median redshift to have higher ratios. 

4.2 Far-infrared SEDs 

Flux boosting of nearly the same magnitude as that experienced 
in the BLAST wavebands ( £12.11 ) also influences our Spitzer mea- 
surements and, to a lesser extent, those made using LAB OCA. 
We therefore apply deboosting factors appropriate for the BLAST 
250-/im sources - appropriate to the level of precision required 
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log 1Q L |R /W 

Figure 6. Td versus deboosted IR luminosity. Galaxies with spectroscopic 
and photometric redshifts are plotted as circles and stars, respectively. A 
radio-loud AGN, selected via its radio morphology, is labelled A'; X-ray 
emitter s are labelled 'X'. T he strong correlation between and L60/Ltm 
seen bv lDunne et al.l fcOOOl) is also evident here. Selection at 250 /im tends 
to find local galaxies with w 20 K, but inevitably also catches more 
luminous, distant galaxies with the same observed and rest-frame values 
of Td that are w (1 + z) X higher. 




Figure 7. (7250 versus redshift, where no K corrections have been ap- 
plied. Circles and stars represent spectroscopic and photometric redshifts, 
respectively. A radio-loud AGN, selected via its radio morphology, is la- 
belled 'A; X-ray emitters are labelled 'X'. The shaded region (around 
z ~ 0.56 ± 0.20, where A res t = 160 /xm for the B LAST 250-/im fil ter) 
represents the ±a range of qieo measured for SINGS foale et al1l2007l) . A 
minor correction has been made to qieo (—0.04) to emulate the effect of 
If -correcting the z = 0.56 radio flux density for the BLAST sample. ±3a 
deviations from the mean value of (7250 are shown by dashed lines. 



(Eales et al. 2009) - to Sir, which we define to be the flux ob- 
served between rest-frame 8 and 1,000 /xm. The large uncertainty in 
the deboosting factor is propagated in quadrature with other uncer- 
tainties. Sir could be determined via modified blackbody fits to the 
70-, 160-, 250-, 350-, 500- and 870-/im flux measurements with (3 
fixed to 1.5 (for example). However, with such well-sampled SEDs 
at our disposal, we concluded that Sir is better determined by in- 
terpolating between the measured flux densities at 24-870 /xm, and 
extrapolating outside the data range to rest-frame 8 and 1,000 /im 
with spectral indices of —1.8 (appropriate for an M82-like SED) 
and +3.5, respect ively, then deboosting according to the recipe in 
lEales et all d2009h . These are the values listed in Table Fig. [5] 
shows the deboosted IR luminosity as a function of redshift for our 
sample. 

We do use modified blackbody fits to estimate the effective 
dust temperature, Td, and these are also presented in Table [2] We 
note that several dust components with a range of temperatures 
would be required to replicate the SEDs accurately over the full 
wavelength range. The average observed-frame Td is 20.5 ± 3.5 K. 
Fig. [6l shows the relationship between rest-frame Td and Lir. 
We expect and see a s trong correlation, fol lowing the work by 
iDunne et all d2QQQh and lChapman et all <2005h : at 250 /xm we are 
sensitive to a particular observed temperature regime; rest-frame 
Td is higher (by a factor 1 + z) than observed Td, so for more dis- 
tant detected sources we expect hotter dust. Lir scales with a high 
power of Td, yielding the strong correlation evident in Fig. [6] 



4.3 The FIR/radio correlation 

We have explored trends in the FIR/radio correlation, both the 
monochromati c version, e.g., £70 = log 10 (ftoMm /Si ,4oomh z ) 
as defined by lAppleton et al.l (120041) - where if -correcting both 
S70/im and Si,4oomh z requires assumptions about SED shape if 
a significant range of redshift is present - and the correlation be- 



tween bolometric Sir and if -corrected radio flux density, where 
both quantities are well constrained. 

We adopt a standard routine so that we can compare like with 
like. First, we exclude a radio-loud AGN, identified on the basis of 
its radio morphology (BLAST J033 130-275604). We do not ex- 
clude the X-ray emitters (see TableQ]). Next, we calculate the error- 
weighted mean and standard deviation of the relevant q parameter, 
noting points that deviate from the mean by more t han 3 a (these 
being 'radio-excess AGN' candidates, as described bv lDonlev et al.l 
120051) . The resultant statistics (and notes) are reported in Table [3] 

4. 3. 1 Monochromatic correlations 

The monochromatic relationship between S250/xm and Si,4oomh z , 
<?250 = log 10 (S 2 50/xm/Si,4ooMHz) is shown in Fig. |7]as a func- 
tion of redshift, and the resultant statistics are reported in Table [3] 
Also shown as a grey box in Fig.|7]is the ±la range of qieo mea- 
sured for the 57 members of the Spitzer IR Nearby Galaxies Survey 
SINGS) sample which are detected at both 160 fim and 1,400 MHz 
Dale et al.ll2007h . At z « 0.6, where the BLAST 250-/im filter 
samples rest-frame 160-/xm emission, values of qieo and #250 are 
entirely compatible with the SINGS average, with or without the 
small radio K correction. We interpret this as evidence for a low 
rate of evolution, at least out to z ~ 0.6. 



4.3.2 K -corrected monochromatic correlations 

Next, we looked at the effect of 'if correcting' the monochromatic 
relationship between 250-/xm and 1,400-MHz flux densities, where 
if is the function that transforms a quantity observed at some fre- 
quency into an equivalent measurement in the rest frame of the 
source being observed. Taking the simple case of a power-law SED, 
where Lff and S° bs are the if -corrected luminosity density and 
observed flux density, respectively, and d\ um is the luminosity dis- 
tance (in m), 
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Figure 8. K corrections for the Spitzer, BLAST and LABOCA wavebands, 
assuming an M 82 SED template, calculated using accurate filter transmis- 
sion profiles and defined such that K = /S° hs , so excluding the (1+z) 
term that takes care of the extra width of the differential frequency element 
in the observed frame (equation 0. 




100 1000 
Rest-frame wavelength Atrn 

Figure 9. Template SED with arbitrary flux scale built from the normalised 
photometry of those 250-/im-selected galaxies with secure identifications 
using the technique described bv lPope etal] (2006). The solid line is a 36- 
K modified blackbody. To create a template from this SED, suitable for 
use in a X-correction algorithm, we interpolated between the means then 
extrapolated to shorter and longer wavelengths with a = —1.8 and +2.25, 
respectively. 



Table 3. Summary of q statistics. 



Type 


Mean 


Standard 


K correction(s) 


of q 




deviation 


or stacking parameters 


250- /j,m- selected 


galaxy sample: 




9250 


2.26 


0.35 


None 


9250 


1.85* 


0.43 


M 82 SED, ol - -0.8 


•3250 


1.82* 


0.43 


Arp220SED, ol = -0.8 


9250 


2.08 


0.34 


Fig.[9]SED, a = -0.8 


9250 


1.94 


0.34 


M 82 SED, measured ol 


9250 


1.91 


0.34 


Arp 220 SED, measured ol 


9250 


2.18 


0.28 


Fig.[9]SED, measured ex. 


91 R 


2.40 


0.29 


None 


91 R 


2.37 


0.31 


a = -0.8 


91 R 


2.41 


0.20 


Measured ol 


SINGS sample (Dale et al. 2007): 




9160 


2.68 


0.41 
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t With M 82-like power-law SED (a - 
* One radio-excess AGN present. 



-1.8) shortward of the peak; 



K converts S° JJoomhz ~ radiation emitted at 2,800 MHz and seen 
at 1,400 MHz - into ^4 00MHz , the 1,400-MHz flux density in the 
rest frame of the emitter. Using the definitions of z and spectral 
index, a, i.e. (1 + z) = z^-est/Vobs and S u oc we find that 

K = S?/S° hs = (l + z)- a ,so 



47rdi 2 um 
(1 + *) 



: 4ttC 



(i + zY+ a ' 



(3) 



For the radio spectrum we use the measured slope between 
610 and 1,400 MHz, adopting a = —0.8 where this is consistent 
with the 610-MHz upper limit, or where no measurement is avail- 
able. This results in a median a of —0.4. In the FIR waveband, 
we explore d two options. O ne is to adopt the observed SED of 
M 82, which Iba ret alJ d2008h found gave the least scatter in #24 . K 
correction^] appropriate for an M 82 SED template in the Spitzer, 
BLAST and LABOCA wavebands, calculated using accurate filter 
transmission profiles, are shown in Fig. [8] We also tried the SED 
of Arp 220 and an SED built from those 250-/im- selected galaxies 
with secure identific ations and r edshif ts, normalised at 250 /im us- 
ing the technique of lPope et al.1 (120061) - see Fig. [9] The latter SED 
made an appreciable difference to the absolute value of ^250 and 
resulted in less scatter. Relevant statistics for ^250 under each set of 
assumptions are reported in Table [3] 



4.3.3 K -corrected bolometric correlations 

Moving from the monochromatic correlation to bolometric quanti- 
ties, where a K co rrection is required only in the radio waveband. 
iHelou et all (Il985l) defined gm such that 



47Tdi 2 um 

(1+*) 



KSZ 



(2) 



where the [1 + z) term takes care of the extra width of the differen- 
tial frequency element in the observed frame. Now imagine we are 
working on a sample of z = 1 galaxies at 1,400 MHz: the factor 



qm - log 



Sfir/3-75 x 10 1 



log 



Si, 



(4) 



Wm- 2 J "VWm^Hz- 1 
where they defined Sfir as the flux between 42.5 and 122.5 pm 

4 Our iDL-based K-correction code, together with SED templates and filter 
profiles, are available on request. 
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lo SlO 44GHz / W HZ ~ 1 

Figure 10. Upper: qm as a function of redshift, using if -corrected radio 
luminosities based on measured radio spectra. The shaded area represents 
±<t; dashed lines are at ±3cr. Lower: Lm versus L J 400MHz • A radio-loud 
AGN, selected via its radio morphology, is labelled 'A ; X-ray emitters are 
labelled 'X'. 



and 3.75 x 10 is the frequency at their mid-band, 80 /mi. Working 
at 42.5-122.5 fim would exploit only one tenth of the wavelength 
coverage now at our disposal so, although we adopt this definition, 
we use Sir - the flux between rest-frame 8 and 1,000 /im - instead 
of Sfir. 

We if -correct the radio flux densities as described in the 
discussion of monochromatic relations. The results are shown in 
Fig. [lOl and reported in Table [3] Utilising all the data at our dis- 
posal, we find qm = 2.41 ± 0.20 for our 250-/mi- selected galaxy 
sample, with no evidence of redshift dependence. 

4.4 Assessing redshift evolution of the FIR/radio correlation 
by stacking 

IMarsden et all b009h . IPascale et al] d2QQ9h and IPatanchon etaD 
(12009 ) have demonstrated the power of stacking analyses based on 
confusion-limited maps and the very steep BLAST source counts. 
We have applied those same techniques to study the evolution of 
#250 and qm by stacking at the positions of the thousands of mid- 
IR- selected galaxies with photometric redshifts in ECDFS. We thus 
eliminate the biases introduced by starting from a flux-boosted and 
confused FIR- selected sample and the problems of demanding se- 
cure identifica t ions i n a densely populated 24-/mi image. Since 
Marsden et al. ( 2009) has shown that the BLAST catalogue only 



comprises «10 per cent of the total flux in the maps, stacking, 
which uses all of the pixels, has the potential to greatly increase 
the sensitivity of the analysis. 

However, some new problems are introduced. We have to 
worry about selection effects due to the catalogue of positions 
upon which we stack - in this case a 24-/xm- selected sample, 
with flux densities determined using IRAC positions as priors 
(Mag nelli et al]2009L with the redshift distribution s hown in Fig. [4] 
and w ith a non- trivial if correction (see Fig.[8}. §5 of lPascale et al.l 
(2009) pr ovides a practical d emonstration of the paradigm de- 
scribed by IMarsden et al.l d2009h and shows that our fundamental 
approach to stacking is accurate to the level of precision required 
here. 

Median- stacking in the radio regime should reduce the conse- 
quen ces of radio-loud AGN be coming more prevalent at high red- 
shift founlop & Peacocklll990l) . Stacking should, therefore, allow 
us to probe the evolution of q for galaxies more representative of 
the general population than is possible for a sample requiring both 
FIR and radio detections, modul o 24-/mi selection b iases. 

We adopt the technique of llvison et "all (l2007ah . dividing the 
Mag nelli et al.l catalogue into six redshift bins, evenly spaced in 
logio (1 + z) between z = and 3, containing 561, 1,441, 2,205, 
1,823, 1,900 and 372 galaxies, respectively. Error- weighted mean 
postage-stamp images are obtained in all the FIR filters, allowing 
us to assess the level of any background emission. Median images 
are made from the VLA and GMRT stacks, to reduce the influ- 
ence of radio-loud AGN. In the radio regime, where the spatial 
resolution is relatively high, making images allows us to conserve 
flux density that would otherwise be lost due to smearing by astro- 
metric uncertainties and finite bandwidth (chromatic aberration) at 
the cost of larger flux density uncertainties. Total flux densities are 
measured using Gaussian fits within AXVS; we determine the ap- 
propriate source centroid and width by making radio stacks using 
every available 24-/xm source (SNR ~ 50), then fix these parame- 
ters to minimise flux density uncertainties and avoid spurious fits; 
at 1,400 MHz, where both forms of smearing are most prevalent, 
flux densities were ~2x higher than the peak values. Such losses 
are not expected where the astrometry is accurate to a small fraction 
of a beam, as it is at 70-870 /xm. 

The top and middle panels of Fig. Qj] show the evolution with 
redshif t of the overall SEP an d the median radio spectral index, a, 
for the iMagnelli et alj (120091) sample. We see that the rest-frame 
8-100-/im portion of the SED becomes progressively more impor- 
tant - in terms of its relative contribution to Sir - as we move to 
higher redshifts, despite the 87 0-^m flux density inc reasing slightly 
with redshift (as expected - iBlain & Longairj|l993l) . We find no 
evidence of significant deviations from the error-weighted mean 
spectral index, o^ 00 = -0.75 ± 0.06, nor of a significant trend 
in a4io 00 across < z < 3. A weighted least-squares fit yields 
a^ 4 00 oc(l+ 2 )°- 14±0 - 20 . 

The lower panel of Fig.fTTIshows the evolution with redshift of 
#250 and q\R, the former if -corrected using the SED template de- 
scribed in ^4.31 and shown in Fig.[9j both if -corrected using mea- 
sured radio spectra. Mean values of both #250 and qm lie within 
1 a of the mean values found for our 250-/mi- selected galaxies (Ta- 
ble[3]). We find, however, that qm evolves with redshift, even across 
< z < 1 where incompleteness should not be a major issue 
in our 24-/xm sample. A weighted least- squares fit suggests qm 
oc (1 + z)~° - 15±0 - 03 . #250 is also found to evolve, though the form 
of the evolution is strongly dependent on the choice of SED adopted 
for if correction, i.e. on the shape of the template SED shortward 
of 200 /xm, emphasising the importance of good spectral coverage. 
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Figure 11. Top: ^ bs^ bs versus rest wavelengt h in six redshift bins , 
meas ured at the p ositions of the 24-/im sample (Magnel li et aD 120091: 
IPascale et"al]|2009l) . Middle: Radio spectral index, a, in the six redshift 
bins. The mean value is —0.75 ±0.06 and there is no evidence of significant 
evolution across < z < 3. Bottom: (7250 and giR in the six redshift bins. 
Dashed line: weighted least-squares fit of the form qiR oc (1 + z) -015 ; 
shaded area: ±1<j prediction of the lLacev et al. | (| 2008) galaxy-formation 
model and the spectrophotometric model of Silva et al. (1998) for galaxies 
with L IR > 10 12 L©, as implemented by Swinbank et al. (2008), updated 
to reflect our wavelength coverage and definition of giR. K-corrections are 
based on the measured radio spectra and - for #250 only - the SED dis- 
cussed in ^4.3l and shown in Fig. [9] The scale is chosen to match Fig. [12] 



The evolution is even stronger for M 82 and Arp 220 SEDs. The 
different behaviour with respect to redshift seen in Figs [TUl and ITTI 
is presumably due to selection biases: the lack of evolution seen in 
Fig. [TO] reflects the FIR sample selection employed there; the evo- 
lution seen in Fig.QTJreflects the IR and radio characteristics of the 
increasingly luminous 24-/xm emitters upon which we are stacking 
at progressively higher redshifts. 

The use of median radio images should exclude radio-loud 
AGN at least as effectively as the 3 -a clip employed in £14.31 so 
the influence of the anticipated evolution (with redshift) of radio- 
loud AGN on the stacked values of q should be minimised. Using 
means instead, which we view as an inferior approach, the aver- 



age spectral index, a 6 



falls to —1.12 as would be expected if 



steep- spectrum AGN contaminate the stacks. This, together with 
the ~3x higher mean radio flux densities, shifts qm by —0.67, al- 
though the form of the redshift evolution remains similar, with qm 
oc(l + ^)- al4±0 - 04 . 

Responding to reports that the radio background is signifi- 
cantly bright er than the cumulative intensity seen i n discrete ra- 
dio em itters (iFixsen et al.ll2009l : ISeiffert et alJl2009b . ISingal et all 
(2009) have speculated that we should see the FIR/radio correla- 
tion evolve, driven by ~0.01-10-/iJy radio activity amongst ordi- 
nary star-forming galaxies. Since this is the flux density regime we 
are probing with our stacking analysis, one might hope to catch 
a glimpse of such evolution. Indeed, the basic f orm of the e volu- 
tion we see in qm is consistent with the idea of ISingal et all and 
our spectral ind ex is consistent wit h that of the radio background 
(a4 ° 00 ~ -0.6,lFixsenetalJl2009h . 

ISwinbank et all d2008h used the GALFORM semi-analytical 
galaxy-formation model to predict mild evolution of qm for lumi- 
nous starbursts due to the ~5-Myr timelag between the onset of star 
formation and the re sulting supernovae, in keeping with the mild 
evolution claimed by Kovacs et al. I d2006h on the basis of 350-/im 
observations of 15 SMGs. Kovacs et all suggested that a change of 



the radio spectral index might account for the deviat ion in q. Revis- 
iting the GALFORM work of ISwinbank et al.l d2008h - calculating 
Sir and qm for galaxies with Lir > 10 12 L in a manner consis- 
tent with this work - results in evolution of the form shown by the 
shaded area in Fig. Qj] (lower panel), i.e. the evolution progresses 
in the opposite sense to that seen, with quantitative agreement only 
in the 2 < z < 3 regime. 

The evolution seen here has implications for all the areas that 
rely on the FIR/radio correlation, e.g. the median redshift of SMGs 
determined via the oft-used Ssso/xm/Si^ooMHz flux density ratio 
would be lower. It will be interesting to discover whether the evo- 
lution we report is confir med by Herschel in the deepest multi- 
frequency radio fields (e.g. lOwen & Morri son 2008). 

Taking the average properties of the 24-/im- and 250-/xm- 
selected samples to be representative, so qm = 2.54 ± 0.26, we 
can suggest a simple transformation between radio flux density 
and SFR, appropriate for samples in which radio-loud AGN con- 
stitute only a small fraction: a radio flux density of Si,4oomh z 
(Wm^Hz- 1 ) can be converted into a bolometric IR flux (rest- 
frame 8-1,000 11m in via: 



a -1 n 12.57+(2.54±0.26) 

OIR = 1U Ol, 400MHz, 

and thence into a star-formation rate (SFR) via: 



(5) 



SERGIO" 10 ^^ Moyr- 1 , (6) 
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stacking into the IR/submm images at the positions of faint radio emit- 
ters, in seven flux density bins, spaced evenly in logio ^i,400MHz • Median 
flux densities in each bin are reported in the key. Bottom: (7250 and qm in 
the seven flux density bins. K corrections were not applied because com- 
plete redshift information was not available. The scale was chosen to match 
Fig.[Hl( bottom). 



where ^ is 1.7 for a Salpeter IMF (dN(m)/dm - 2.35) cov 



ering 0. 1-100 Mp 
Scovill 
2000h . 



10 (Condon 1992; Kennicut 
Scoville & Young l 1 19831 : iThronson & Telesc 



1998al l 3, see also 
19861 : Tlnoue et all 



4.5 Stacking at the positions of sub-mjy radio galaxies 

The value of qm found in the previous section is weighted by the 
types of galaxies found in the 24-/im catalogue. To investigate 
whether qm is different for a radio- selected population, we stack 
the FIR/submm images at the positions of 816 robust radio emit- 
ters (>5 a, 40 < Si, 4oomhz < 1, 000 /iJy), divided into seven flux 
density bins, spaced evenly in logio *Si,4oomh z . 

The FIR SEDs of the radio- selected galaxies are qualitatively 
similar to those of low-redshift 24-/xm- selected galaxies (compare 
the top panels of Figs Qj] and [12]). We find that ^250 and qm are 
lower for radio- selected galaxies than for those selected in the mid- 
IR or FIR wavebands (Table [3]), as one might expect given the se- 
lection criterion. However, the difference is surprisingly small: for 
40 < Si, 4oomh z < 100 pjy we obtain q 2 5o = 2.42 and qm = 2.25, 
cf. the equivalent (pre-if -correction) values for our FIR-selected 
galaxy sample, ^250 = 2.26 =b 0.35 and qm = 2.40 =b 0.29. Bear 



in mind that, at a redshift of unity, qm would rise by ~0.06 due to 
the radio K correction and fall by ~0.15 due to the shift of A res t 
= 8-1000 pm to lower observed frequencies - an overall shift of 
~ —0.1. This suggests that star formation plays a significant role 
in powering faint radio galaxies, providing over half of their IR 
luminosity. Radio selection will inevitably have favoured galaxies 
with recently injected relativistic electrons so this is likely to be a 
lower limit. We find that ^250 and qm anti-correlate with radio flux 
density (Fig. Q21 lower panel) and we suggest this is due to the in- 
creasing pre valence of radio-loud AGN as one moves into the mJy 
regime (e.g. llbar et ap2009h . 



5 CONCLUSIONS 

We have defined a sample of 25 0-/xm- selected galaxies using data 
from BLAST. The noise is dominated by confusion, which severely 
limits the number of robust detections and hampers the identifica- 
tion of secure, unambiguous mid-IR or radio counterparts. These 
problems were not always apparent when using an objective, prob- 
abilistic approach to cross-identifying the galaxies. 

We find that the most likely 24-/xm counterparts to our small 
sample of 250-/xm galaxies have a median [interquartile] redshift 
of 0.74 [0.25, 1.57]. 

At z w 0.6, where the BLAST 250-/im filter probes rest- 
frame 160-/xm emission, we find no evidence for evolution of (/250 
relative to qiQo measured for the SINGS sample of local galaxies. 

We find that Sir is better correlated with radio flux density 
than S250/xm, and that if -correcting the radio luminosity using 
measured spectral slopes reduces the scatter in the correlation. 

qm - the logarithmic ratio of bolometric IR and monochro- 
matic radio fluxes - is determined for FIR-, mid-IR and radio- 
selected galaxies. We provide a simple recipe to convert a radio 
flux density into an instantaneous SFR, via bolometric IR luminos- 
ity, for mid-IR- or FIR-selected samples that do not contain a large 
fraction of radio-loud AGN. 

qm is found to be similar for our 250-/xm- and radio- selected 
galaxies, which suggests that star formation is responsible for over 
half of the IR luminosity in the latter, especially the faintest ra- 
dio galaxies (Si,4oomh z < 100 /iJy). This fraction could well be 
higher given that radio selection favours galaxies with recent injec- 
tions of relativistic electrons. 

Stacking into 610- and 1,400-MHz images at the positions 
of 24-/im- selected galaxies, we find no evidence that the spectral 
slope at radio wavelengths evolves significantly across < z < 3. 
However, we find tentative evidence that qm does evolve, even 
across < z < 1 where incompleteness in the parent sample 
should not be a serious issue. The evolution is of the form, qm 



oc (1 + z) 



, across the peak epoch of galaxy formation. 



This has major implications for many techniques that rely on the 
FIR/radio correlation. We compare with semi- analytical model pre- 
dictions and speculate that we may be seeing an increase in radio 
activity amongst ordinary, star-forming galaxies - relative to their 
IR emission - amon gst that has been suggested may give rise to the 
radio background dFixsen et al.ll2009l : ISingal et alJ2 009). 

It will remain difficult to associate galaxies detected in Her- 
schellSPlKE surveys unambiguously with IR or radio counterparts, 
since the modest increase in telescope aperture and the steep source 
count s will ensure that confusion remains an issue toevlin et al.1 
2009). However, alongside 70-, 100- and 160-/im imaging from 
Herschel/PACS and 450- and 850-/zm SCUBA-2 data from the 15- 
m James Clerk Maxwell Telescope, it should be possible to fine- 
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tune the astrometry to a ~ 1 arcsec and to deblend sources in deep 
SPIRE images. For the ~ 1-2 deg 2 of FIR/submm coverage planned 
for the fields with the deepest available radio and Spitzer data, e.g. 
GOODS, the Subaru/XMM-Newton Deep Survey and the Lock- 
man Hole, with ai,4ooMHz ^ Soybean! -1 and cr24um ^ 10/iJy 
(iBiggs & IvisodbOOd lOwen & Morrisonll2008l : iMiller et al Eooil : 
Morrison et al., in preparation; Arumugam et al., in preparation), 
almost every deblended FIR source should have relatively secure 
24-/xm and radio identifications, the majority unambiguous, pro- 
viding a complete redshift distribution and less biased estimates of 
qm for FIR- selected galaxies. 
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